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We analyze electronic and phononic quantum transport through zigzag or chiral graphene nanorib- 
bons (GNRs) perforated with an array of nanopores. Since local charge current profiles in these 
GNRs are peaked around their edges, drilling nanopores in their interior does not affect such edge 
charge currents while drastically reducing heat current carried by phonons in sufficiently long wires. 
The combination of these two effects can yield highly efficient thermoelectric devices with maximum 
ZT ~ if at liquid nitrogen temperature and ZT ~ 4 at room temperature achieved in ~ f (im long 
zigzag GNRs with nanopores of variable diameter and spacing between them. Our analysis is based 
on the 7r-orbital tight-binding Hamiltonian with up to third nearest-neighbor hopping for electronic 
subsystem, the empirical fourth-nearest-neighbor model for phononic subsystem, and nonequilib- 
rium Green function formalism to study quantum transport in both of these models. 
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The recent explosion of research on graphene — one- 
atom-thick allotrope of carbon — has been largely focused 
on its unique electronic structure and transport proper- 
ties governed by the two-dimensional honeycomb lattice 
of carbon atomsP Very recently, the exploration of its 
thermal and thermoelectric properties has been initiated 
by measuring the thermopowei'2lEl S and phonon thermal 
conductivity^ if ph of large-area graphene. The measured 
values^ of S ~ 100 juV/K near the Dirac point (DP), 
as well as the room-temperature K p h ~ 4000 W/mK 
(averaged over values obtained using different samples 
and experimental techniques^) which outperforms virtu- 
ally all other known materials, point out that large-area 
graphene is not suitable for thermoelectric applications. 

Thermoelectrics transform temperature gradients into 
electric voltage and vice versa. Although a plethora of 
widespread applications has been envisioned, their usage 
is presently limited by their small efficiency.'- Thus, care- 
ful tradeoffs are required to optimize the dimensionless 
figure of merit 



which quantifies the maximum efficiency of a thermoelec- 
tric cycle conversion in the linear-response regime where 
small voltage V = —SAT exactly cancels the current 
induced by the small thermal bias AT. This is due to 
the fact that ZT contains unfavorable combination of 
S, average temperature T, electronic conductance G and 
thermal conductance n c \ + /c p h. The total thermal con- 
ductance has contributions from both electrons n e \ and 
phonons K p h- The devices with ZT > 1 are regarded as 
good thermoelectrics, but values of ZT > 3 are required 
for thermoelectric devices to compete in efficiency with 
conventional power generators and refrigerators P 

Thus, a number of proposals have been put forth to 
evade the problem of high lattice thermal conductivity 
of large-area graphene that could open a pathway for 
its thermoelectric applications. For example, large-area 



graphene could reach ZT s» 0.3 if perforated by the so- 
called antidot lattice tailored to impede phonon prop- 
agation. 7 Switching to quasi-one-dimensional graphene 
nanoribbons (GNRs) makes possible further enhance- 
ment of ZT where it has been predicted that long (~ 1 
/im) GNRs with zigzag edges and disorder introduced 
along such edges could reach ZT ~ 4 at room tempera- 
tureP Another route is to engineer structural defects in 
GNRs that can block phonons while retaining quasibal- 
listic electronic transport!^ 

However, it is more advantageous to search for high- 
ZT devices among clean nanowired^ since edge or sur- 
face disorder can affect electronic conductance signifi- 
cantly. For example, the experiments on etched GNRs 
with rough edges find Coulomb blockade effects (not 
taken into account in Ref. |5]) and transport gap much 
larger than the band gapFEQ 

In this Rapid Communication, we exploit peculiar elec- 
tronic transport properties of clean GNRs with zigzag 
(ZGNR) or chiral (CGNR) edges, illustrated in Fig. [TJ 
where the local charge current density carried by quasi- 
particles sufficiently close to the DP is peaked around 
nanoribbons edges as demonstrated in Fig. [2] Thus, 
drilling nanopore*^ in the ZGNR or CGNR interior will 
not substantially modify such "edge currents" . This is 
confirmed by the transmission function in Figs. [3]ja) and 
j3^c) which is reduced from T e i{E) = 3 in homogeneous 
GNRs to %i(E) ~ 2 around the DP for both ZGNR 
and CGNR with an array of nanopores. Furthermore, 
Tei(E) around the DP does not change as one increases 
the length of GNRs because "edge currents" propagate 
quasiballistically. 

On the other hand, nanopore arrays, whose fabrica- 
tion has been pursued recently by a variety of experi- 
mental techniques,^ break homogeneity of the nanowire 
so that they can substantially impede the propagation of 
phonons in sufficiently long GNRs. This is corroborated 
by our results for the phonon transmission function in 
Figs. Qa) and ^b) and the corresponding lattice ther- 
mal conductance in Figs.^c) and Etcl). 
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FIG. 1: (Color online) Schematic view of: (a) 20-ZGNR (com- 
posed of 20 zigzag chains); and (b) (8,1)-CGNR with chiral 
angle = 5.8°. The size of the nanopores, assumed to be 
drilled in the GNR interior away from its zigzag or chiral 
edges, and the distance between them is illustrated by plot- 
ting two repeated supercells of each GNR. The length of these 
GNRs in actual calculations is set to L ~ 1.2 /im. 



Combining these two effects, we obtain maximum 
ZT ~ 3 at T = 77 K and ZT ~ 1.5 at T = 300 K 
in Fig. |5jc) for the case of (8,1)-CGNR in Fig. [I] whose 
identical nanopores are arranged in a perfectly symmet- 
ric and ordered array. The values of ZT for 20-ZGNR 
with periodic array of nanopores are lower, as shown in 
Fig.J^a]. In realistic GNR-based devices, it may be chal- 
lenging^ to control the pore arrangement to a high de- 
gree of order shown in Fig. [T] If we assume that pore 
diameter takes a random value within some interval and 
that nanopores are not arranged into perfectly periodic 
array, then we find in Fig. 5jb) a possibility of enormous 
ZT ~ 11 at T = 77 K and ZT ~ 4 at T = 300 K in the 
case of ZGNR-based device. 

In the rest of the paper we explain details of our mod- 
els for electronic and phononic subsystem and the corre- 
sponding quantum transport calculations. The early the- 
oretical studies of ZGNR-based devices have utilized^ 
a simplistic tight-binding model (TBM) with single tt- 
orbital per site and the nearest-neighbor hopping only, 
or its long-wavelength (continuum) approximation — the 
Dirac-Weyl Hamiltoniari™ — valid close to the DP. How- 
ever, both of these models predidPE^ that transmission 
function of an infinite homogeneous ZGNR is T e \ = \ 
around the DP and that current density profile is 
peakecf^in the middle of ZGNR (even though local den- 
sity of states reaches maximum around the edges). This 
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FIG. 2: (Color online) Spatial profile of local charge cur- 
rents over the transverse cross section of infinite homoge- 
neous 20-ZGNR or (8,1)-CGNR for electronic transport close 
(E F = -0.43 eV) to the DP. The sum of bond current^ 1 
Jnm/V, which describe charge flow from site n to site m of 
the honeycomb lattice if hopping t™ 7^ is non-zero between 
the two sites, gives the conductance G = I/V (I is the total 
current in the leads and V — > is small bias voltage driving 
the linear-response transport). 



contradicts first-principles calculations^, or TBM with 
up to third nearest-neighbor^ hopping parameters fitted 
to such first-principles calculations, which predict T c \ = 3 
around the DP, as well as that local current density is 
mostly confined to flow around the zigzag edges ff^ It is 
worth mentioning that the majority of recent studies of 
thermoelectric properties of ZGNRs with edge disorder^ 
or of finite length graphene antidot latticed have utilized 
the TBM with nearest-neighbor hopping, so that a pos- 
sibility to exploit "edge currents" around zigzag or chi- 
ral edges for thermoelectric device applications has been 
overlooked. 

Most importantly, the recent experiments have con- 
firmed the existence of "edge currents" in metallic 
ZGNRs by actually utilizing them to increase the heat 
dissipation around edge defects and, thereby, rearrange 
atomic structure locally until sharply defined zigzag 
edge is achieved PS Also, the very recent chemical syn- 
thesis^! of (8,l)-CGNRs via carbon nanotube unzip- 
ping method have exhibited properties in sub-nanomctcr- 
resolved scanning tunneling microscopy and spectroscopy 
that can only be explained by the existence of smooth 
edges supporting edge quantum states (i.e., wavefunc- 
tions whose probability density is l arg e around the 
edges). Although ZGNRs or CGNRspTare insulating 
at very low temperatures due to one-dimensional spin- 
polarized edge states coupled across the width of the 
nanoribbon, such unusual magnetic orderi ng an d the 
corresponding band gap is easily destroye d 21 * 22 * above 
T > 10 K. Therefore, both ZGNRs and CGNRs can be 
considered as metallic nanowires at liquid nitrogen or 
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FIG. 3: (Color online) (a) Zero-bias electronic transmission 
%\{E) for an infinite homogeneous 20-ZGNR or 20-ZGNR of 
length L ~ 1.2 /im with periodic array of identical nanopores 
shown in Fig. [T^a). (b) Thermopower at two different tem- 
peratures for finite length 20-ZGNR with nanopores. (c) 
Zero-bias electronic transmission for an infinite homogeneous 
(8,1)-CGNR or (8,1)-CGNR of length L ~ 1.2 with pe- 
riodic array of identical nanopores shown in Fig. |TJb) . (d) 
Thermopower at two different temperatures for finite length 
(8,1)-CGNR with nanopores. 

room temperature analyzed in our study. 

We adopt the TBM with single 7r-orbital per site and 
up to third nearest-neighbor hopping 

H = J2 £ " fi n £ " - E ^Pm. (2) 
1 1 n,m 

to describe the electronic subsystem of 20-ZGNR and 
(8,1)-CGNR in Fig.[l] The operators (c n ) create (an- 
nihilate) electron in the 7r-orbital located on site n of the 
honeycomb lattice whose lattice constant is a s» 0.246 nm 
and C-C bond length is d^g w 0.142 nm. For clean 
GNRs studied here the on-site potential is set to zero, 
e n = 0. We consider up to third nearest-neighbor^ hop- 
ping parameters— t^ +dAB = 2.7 eV, t^ +dAA = t^ +dBB = 
0.2 eV, and i£ +dAB ' = 0.18 eV— which describe the 
nearest-, next-nearest- and next-next-nearest neighbor 
hopping, respectively. Since the honeycomb lattice of 
graphene is composed of two triangular sublattices A and 
B, the parameters t" +dAB and i£ +dAS ' describe intersub- 
lattice hopping, while t™ +dAA = tJ +dBB describes the 
intrasublattice hopping. 

In the coherent transport regime, the retarded Green 
function (GF) matrix^ 

G(E) = [EI-H-V L (E)-X R (E)}-\ (3) 

makes it possible to expresses the zero-bias electron 
transmission function between the left (L) and the right 



FIG. 4: (Color online) (a) The phonon transmission func- 
tion 7ph(w) and (b) the corresponding phonon thermal con- 
ductance K ph for an infinite 20-ZGNR or 20-ZGNR of length 
L ~ 1.2 pirn with periodic array of identical nanopores shown 
in Fig. [TJa) . (c) The phonon transmission function and (d) 
the corresponding phonon thermal conductance K p h for an in- 
finite (8,1)-CGNR or (8,1)-ZGNR of length L ~ 1.2 (jm with 
a periodic array of identical nanopores shown in Fig. |TJb) . 

(R) electrodes as: 

T ei (E) = Tr {T R (E)G(E)T L (E)G^(E)} . (4) 

Here H is the matrix representation of Hamiltonian in 
Eq. ([2]) , I is the unit matrix in the Hilbert space of the 
active region, T,l :R (E) are the self-energies due to the 
"interaction" with the leads, and Tl v h(E) = i[£L :R (E) — 
H\ R {E)\ are the level broadening matrices determining 
the escape rates for electrons to enter into the attached 
leads. In realistic devices, active region consisting of 
ZGNR or CGNR of finite length with nanopores will 
eventually need to be connected to metallic electrodes. 
However, since GNR+nanopores devices we analyze are 
rather long ~ 1 /jm, and screening takes place over a dis- 
tance much shorter than the active region, it is justified 
to use semi-infinite homogeneous ZGNRs or CGNRs as 
leads for simplicity. 

The transmission function Eq. Q allows us to compute 
the following integrals^ 

00 

K n Qj) = l J dET c m{E-v) n [- d -^§f^j, (5) 

— OO 

where f(E, fi) = {1 + cxp[{E - fij/keT}}- 1 is the 
Fermi-Dirac distribution function at the chemi- 
cal potential /i. The knowledge of K n (fi) finally 
yields all electronic quantities in the expression for 
ZT: G = e 2 K (fi); S = K^/leT^fj,)]; and 
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FIG. 5: (Color online) The thermoelectric figure of merit ZT 
for GNR+nanopores devices of length L ~ 1.2 /jm: (a) ZT 
vs. energy at two different temperatures for 20-ZGNR with a 
periodic array of identical nanopores illustrated in Fig. [IJa); 
(b) ZT vs. energy at two different temperatures for 20-ZGNR 
with pores of variable diameter D £ [4.5 cLab, 7.5 (Iab\ and 
their position shifted by a value Ax € [— 2dAB,2dAfl] along 
the transport direction away from the original position within 
periodic array shown in Fig. [Tf a) ; (c) ZT vs. energy at two 
different temperatures for (8,1)-CGNR with a periodic array 
of identical nanopores illustrated in Fig.JTJb); and (d) ZT vs. 
energy at two different temperatures for (8,1)-CGNR with 
pores of variable diameter D £ [1.5 cLab, 3.5 cLab] and their 
position shifted by a value Aa; 6 [—0.5 cLab, 0.5 <Iab] along 
the transport direction away from the original position within 
periodic array shown in Fig. |TJb) . 



The phonon subsystem is modeled using the empirical 
fourth-nearest-neighbor model. Its parameters were orig- 
inally fitted to match the phonon dispersion measured by 
Raman spectroscopy^ and X-ray scattering data. Our 
model is reparametrizecP^ to include the off-diagonal 
terms of the force constant matrices and rotational in- 
variance conditions which require a small correction to 
the in- and out-of-plane tangential force. Such empirical 
interatomic potential model offers excellent fit to both 
experiments and first-principles numerics.^ 

The phonon thermal conductance, in the absence of 
phonon-phonorP-O or electron-phonorP^ scattering, is ob- 
tained from the phonon transmission function T p h(uj) us- 
ing the Landauer-type formula^ 



transport regime can be expressed in complete analogy 
with Eq. Q for elastic electronic transport 

TphM = Tr {A^HDHAiO^DV)} . (7) 
The phonon retarded GF is obtained in the same fashion 
as the electronic one in Eq. ([3| but with substitutions 
H -> K, EI -> uj 2 M and -> U r 



D(w) = [uj 2 M K - U L (u) - n R (u)Y 



(8) 



Here K is the force constant matrix, M is a diagonal ma- 
trix with the atomic masses, 11^^ are the self-energies, 
and A L _ R ((jj) = i[TL LiR (uj) — 11^ R {w)]. This methodol- 
ogy does not take into account resistive umklapp phonon- 
phonon scattering which plays an important role in in- 
terpretation of experiments on room-temperature lattice 
thermal conductivity of large-area grapheneP However, 
this effect, which is easy to describe using the Boltzmann 
equation but is very expensive computationally within 
the noncquilibrium GF formalism,^ should not play an 
important role in nanoribbons depicted in Fig. [T] because 
their width is much smaller than the mean-free path 
I ~ 677 nm due to phonon-phonon scattering in large- 
area graphene at room temperature!^! 

In conclusion, quantum transport analysis of electronic 
and phononic transport in ZGNR and CGNR, where 
electronic subsystem is described by the 7r-orbital tight- 
binding Hamiltonian with up to third nearest-neighbor 
hopping and phononic subsystem is described by the 
empirical fourth-nearest-neighbor model, suggests that 
these nanowires could serve as building blocks of highly 
efficient thermoelectric devices when perforated by an 
array of nanopores residing in their interior. This is due 
to the fact that local charge current density is peaked 
around their edges, as demonstrated explicitly by Fig. [2] 
and confirmed experimentally,^ so that nanopores do 
not impede such "edge currents" while drastically re- 
ducing phonon conduction in sufficiently long ZGNRs 
or CGNRs. In the case of periodic array of identical 
nanopores, we find that largest ZT ~ 3 at T = 77 K 
and ZT ~ 1.5 at T = 300 K can be reached using (8,1)- 
CGNR. On the other hand, if the pore diameter takes 
a random value within some interval and the distance 
between the pores is varied, then we find a possibility 
to further optimize the figure of merit which can reach 
astonishingly large values, ZT ~ 11 at T = 77 K and 
ZT ~ 4 at T = 300 K, in the case of ZGNR-based de- 
vice. 
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